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We explore possible signatures for observing Majorana fermions in the tunneling spectroscopy 
of high-T c superconductors. It is shown that due to the Rashba spin-orbit interaction (or) gener¬ 
ally introduced by contact electrode, in addition to the Heisenberg spin exchange interaction, the 
Dzyaloshinskii-Moriya and spin dipole-dipole interactions are induced. As a result, p- wave super¬ 
conductivity is induced with the gap function d- vector being not aligned with the internal magnetic 
field of the spin-orbit interaction. For typical strength of the Rashba interaction, the induced p- wave 
is weak. Hence the resulting superconductor is still gapless and is not a topological superconductor. 

However, we find that the ground state undergoes a phase transition to a topological gapless phase 
with each nodal point originated from pure d- wave being split into two stable nodal points char¬ 
acterized by the symmetry class DHL Due to the splitting nodal structure, zero-energy Majorana 
modes always exist for any interfaces that are not exactly in (100) or (010) directions. Hence for 
general interfaces, existence of Majorana modes is a robust feature. In addition, due to the non- 
aligned d- vector, for (110) interfaces in which d- wave is subdominant to p-wave, there exist sizable 
dispersive Majorana edge states. Our results indicate that due to the presence of these Majorana 
modes, a small plateau in tunneling spectrum near zero bias peak would be induced. Furthermore, 
zero-energy Majorana modes result in 47r periodicity in typical SIS 7 junctions with difference in 
orientations of S and S' being within 21° — 39° for aR — 0.05 eV—0.3 eV. As a result, it is easy for 
a 7r-ring in tricrystal experiments to hold Majorana fermions and exhibit periods of two flux quanta 
in external magnetic field. These phenomena may have been already observed in experiments and 
their connections to experimental results are discussed. 


I. INTRODUCTION 

In the past decades, searching for the mysterious par¬ 
ticle, Majorana fermion, proposed by Ettore Majorana 
in 1937 2 has continued to be an important issue in high- 
energy physics. Recently, the field of topological super¬ 
conductors and superfluids 1 in condensed matter physics 
has been drawing much attention since these topological 
matters are capable of hosting Majorana fermions. In¬ 
stead of seeking for it in high-energy experiments, topo¬ 
logical matters in condensed matter systems may shed 
light on finding this elusive Majorana fermion. Majo¬ 
rana fermion, half state of a fermion, obeys non-Abelian 
statistics and possesses the non-local property. Both of 
the significant ingredients render such topological mat¬ 
ters a great potential to realize quantum computers — 
the dream which scientists have been pursuing for over 
30 years. 

Topological matters are classified by distinct topo¬ 
logical numbers based on dimension and symmetries of 
the system. In analogous to the framework of topo¬ 
logical insulators characterized by a Z 2 invariance, non- 
centrosymmetric superconductors which break inversion 
symmetry but preserve time reversal symmetry are gen¬ 
erally termed as time-reversal-invariant (TRI) supercon¬ 
ductors. Inversion symmetry breaking leads to the mix¬ 
ing of spin-singlet and spin-triplet Cooper pairing^ 3 . As 


a result, intensive studies have been dedicated to investi¬ 
gate fully gapped TRI superconductors for their poten¬ 
tial to be a topological nontrivial matter, termed as TRI 
topological superconductors 4 K In order to have fully 
gapped superconductors, many proposals for creating 
Majorana fermions rely on the proximity effect of 8-wave 
superconductor^. The disadvantage of these approaches 
lies in the fact that most 8-wave superconductors have 
low critical temperatures and Majorana fermions can be 
generated only in very low temperatures. It is therefore 
desirable to seek alternative approaches based on high-T c 
cuprate superconductors. 

In real materials, the heavy-fermion superconductor, 
CePt 3 SP@, discovered in the early days and the recent 
finding of superconductivity at the interface between 
LaAlOs and SrTi 03 10 ^ are typical examples of TRI su¬ 
perconductors in which spin-orbit interaction comes into 
play due to lack of inversion center. Recently, it is ar¬ 
gued that the contact of heavy metals with supercon¬ 
ductors can induce the Rashba spin-orbit interaction in 
the interface 11 . In particular, for Au/YBCO interface, 
it is estimated that the Rashba spin splitting energy at 
Fermi surface is of the order 200 meV“ There is also ex¬ 
perimental evidence showing that there may exist intrin¬ 
sic spin-orbit interaction in cuprate superconductors 13 . 
From theoretical point of view, however, the strength of 
available spin orbit interaction is too small to turn a d- 
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wave cuprate superconductor into a fully gapped TRI su¬ 
perconductor. Therefore, even though the ideal scheme 
is to generate Majarona fermions via TRI topological su¬ 
perconductors, its feasibility is low in high T c cuprates 
and it is necessary to seek for alternative approaches. 


In this paper, we investigate effects of the Rashba 
spin-orbit interaction on high-T c cuprate superconduc¬ 
tors. Based on the mean-field theory of t-J model^H^ 
we show that the gap function d-vector for the induced 
spin-triplet pairing is generally not aligned with the inter¬ 
nal magnetic field of the spin-orbit interaction. This fact 
respects the C± v point group symmetry and for typical 
strength of the Rashba interaction, the superconductor 
is still gapless with each nodal point of the original d- 
wave superconducting state being split. The resulting 
nodal points are protected by symmetries in the symme¬ 
try class DIII^. Furthermore, due to the splitting nodal 
structures, zero-energy Major ana edge modes that are 
determined by projection of bulk nodes to edges 17 always 
exist for any orientations of interfaces that are not exactly 
in [100] or [010] directions. In particular, near [110] di¬ 
rection in which d-wave is subdominant to p-wave, there 
exist sizable dispersive Majorana edge states. Our re¬ 
sults indicate that the presence of these Majorana modes 
gives rises to a small plateau in tunneling spectrum near 
zero bias peak. Furthermore, these Majorana modes al¬ 
ways result in 47r-Josephson effect in typical junctions 
and hence flux trapped in tricrystal experiments jump 
in unit of two flux quanta. These phenomena may have 
been already observed in experiments and their connec¬ 
tions to other experiments will be discussed. 


This paper is organized as follows. In Sec. [TT| by in¬ 
cluding the Rashba spin-orbit interaction, we derive su¬ 
perconducting phase diagram based on the mean-field 
theory of the effective low-energy Hamiltonian for the 


strong coupling limit of the Hubbard model. In Sec. Ill 


the gapless regime of d + p wave resulted from typi¬ 
cal strength of Rashba interaction is investigated. Edge 
states and their relations to nodal point structures are 
examined. In particular, how the edge states of (110) 
edge change in the presence of the Rashba spin-orbit in¬ 
teraction changes is also examined. In Sec. IV, we ex¬ 


amine the tunneling spectrum for the NIS junction and 
periodicity of Josephson current in a SIS' junction. The 
contribution due to the dispersive and flat-band Majo¬ 
rana modes are examined. Furthermore, the relation of 
47r periodicity and flat-band Majorana modes are derived 
for Josephson junctions. The conditions for the 7r-ring in 
the tricrystal experiment to hold Majorana fermions are 
also examined in details. Finally, in Sec. [V| we conclude 
and discuss possible connections of our results to exper¬ 
imental observations. 


edge contact 



FIG. 1: Typical configuration of high-T c cuprate supercon¬ 
ductors in contact with other materials (electrode or sub¬ 
strate). Through the interface between the cuprate and the 
underlying material, spin-orbit interactions that break inver¬ 
sion symmetry are induced in the cuprate superconductor 
either through proximity effect or in the layer near the in¬ 
terface. Here the cuprate superconductor can make further 
tunnel junction with other superconductors via the edge con¬ 
tact. 


II. THEORETICAL MODEL AND MEAN-FIELD 
PHASE DIAGRAM 


The parent compound of high-T c cuprate superconduc¬ 
tor is well-known to be an antiferromagnetic Mott insu¬ 
lator. In spite of the fact that a complete theory for high- 
Tc superconductors is under debate, a common consensus 
for grasping the essential physics behind the emergence of 
superconductivity lies in the 2D Hubbard model within 
large- U limit, or the equivalent t-J modeP^l However, 
in typical experiments, high-T c cuprate superconductors 
are often in the form of thin films with contact to other 
materials along c-axis as illustrated in Fig. [l] Through 
the interface between the cuprate and the underlying ma¬ 
terial, spin-orbit interactions that break inversion sym¬ 
metry may be induced either in the cuprate supercon¬ 
ductor by proximity effect or induced in the layer near 
the interface. In particular, for Au/YBCO interface, it is 
estimated that the energy splitting due to the spin-orbit 
interaction is of the order 0.2 eV 12 . The largest splitting 
reported is also around 0.3e\Eal. It is therefore necessary 
to include effects due to the spin-orbit interaction up to 
the spitting energy being around 0.3eV. 

In the presence of the Rashba spin-orbit interaction, 
however, there is no full understanding of correspond¬ 
ing superconducting phases. From the symmetry point 
of view, the Rashba spin-orbit interaction breaks the in¬ 
version symmetry. As a result, parity is not conserved 
and hence it is generally true that singlet pairing will 
be mixed with the triplet pairing in any superconducting 
phases. To be explicit, we shall start from the 2D Hub¬ 
bard model with the Rashba spin-orbit interaction. The 
effect of dispersion along c-axis will be examined at the 
end of Sec. IV. In the strong coupling limit, exchange in¬ 
teractions give arise an extended t — J model which can 
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be generally expressed as H e ff = H t + Hj with 

Ht = — ^4 ct Cj, a 
ij 

+ ^2 [ CT X ( r i “ ri )la/3 ' ^IcA* 0 ) + 

<ij> 

(1) 

Hj = H Jl + Hj 2 + i^j 3 , 

= ^ Ji MS* • S, - j + J 2 [D • (S* x S^)], 

<L/> ' ' 

S* • Sj + ^riiTij — 2(D • Sj)(D • Sj) . 

( 2 ) 

Here a is the lattice constant. Uj = t, £', and t" for 
sites i and j are nearest, next nearest, and the third 
nearest neighbors, respectively, and Uj = 0 for longer 
distance, ctR characterizes the energy strength of the 
Rashba spin-orbit interaction. Q ?cr = (1 — sat¬ 

isfies the no-double-occupancy constraint. For electrons 
In terms of the on-site Coulomb interaction (Hubbard I/), 

three spin interactions are given by J% = , J 2 = — 

2 

and J 3 = — fhf-, representing the Heisenberg interaction, 
the Dzyaloshinskii-Moriya (DM) interaction and the spin 
dipole-dipole interaction respectively. The orbital vector 
D that couples to spins is related to the lattice constant a 
connecting sites in nearest neighbors along x(y ) direction 
by D = ((rj -r i) y ,-(rj -r^O) = (a,-a, 0). 

To satisfy the no-double-occupancy constraint, one can 
resort to the slave-boson method 14 or adopt Gutzwiller 
approximations by using renormalized parameters 19 . In 
low doping regime, both approaches yield the same mean- 
field Hamiltonian. We shall adopt the Gutzwiller ap¬ 
proximation. In the Gutzwiller approximations, strength 
of spin interactions remain the same in the low doping 
regime. Furthermore, is replaced by c* jCT . The spin 
interactions can be generally factorized as summations 
of particle-particle channels and particle-hole channels 20 . 
To investigate superconducting phases, one focuses on 
particle-particle channels. We find that the J\ term 
yields the spin-singlet pairing channel. By contrast, the 
J 3 term contains the spin-triplet pairing channel, and 
the J 2 term mixes the spin-singlet and spin-triplet pair¬ 
ing channels. Specifically, the pairing channels can be 
written as 

H Ji = ( 3 ) 

Hj 2 = J 2 ]T (D*J • 4-)4 + $,-(Dy • hi), (4) 
<ij> 

Hj 3 = (5) 

<b'> 

where the singlet pairing field ijjji = Cj^c^ — Cj^c^, 
the spin-triplet vector field = (cj^ct^, 0) and 


+ J 3 Y, 


TABLE I: The representations of Ca v 



Even parity 

Odd parity 

Ai 

^4i 9 1 

A 2u k x y - k y x ; 

k x ky(k x x k y y)\ 

k x k y k z (kl - ky)z 

A 2 

A 2g k x k y (kl - kl) 

Aiu k z z ; k x x + k y y 

B 1 

Big kl - kl 

B 2u k x y + k y x 

b 2 

B 2 g k x ky 

B\ u k x x kyy 

E 

Eg k z (k x zb iky ) 

E u ±iz(kl +n ±ikl +n )- 

±ik z (k™x zb ikyy); 
±ik x kyk z (kxii =b ikyX ), 
n = 0,2 


$ 


T = ^ (r>' r>' A s*' Oi 

The vector D.y are given by D^ +x = — | (1,1,0) and 
i | (—1,1,0). In addition to superconducting 
channels, we also include particle-hole chanels and set 

XOl = E A C l C i+aia) and Xll = {4-\ C i+ah) with 1 = X 
or y. When the hole concentration is S, the resulting 
mean-field Hamiltonian Hmf can be written as 


Hmf = ^2 & C Lck* C ka ' Sk) a(r > C ka' 

kcr kcr cr' 

+ (4 fC t_ H - cld kt ) 

k 

+ ( 6 ) 

k 


where the lattice constant a is set to be 1, ^ = 
—2[t(cosk x + cos k y ) + 2St' cos k x cos k y + St" (cos2k x + 
cos2 k y )] — y and = (smk yi — sink x ) with t = St + 

(Jl 2 J3) Xo - if Xix and a R = 6a R - (J x + J 3 )xix ~ %Xo■ 
The pairing amplitudes A s and A t represent amplitudes 
for singlet and triplet pairing with corresponding gap 
functions Y k and Y k respectively. The allowed represen¬ 
tations for gap functions are tabulated in the table [H| 
In the absence of the Rashba spin-orbit interaction, it is 
known that d-wave in B\ g representation and s-wave in 
A\ g representation are two superconducting states with 
lowest energies. In the presence of the Rashba spin-orbit 
interaction, we find that A i and B\ representations be¬ 
come the lowest two in energy due to mixing of the singlet 
and triplet pairings. For s + p pairing symmetries in A\ 
representation, we have Y 9 k = cos k x + cos k y and Y k = 
sin k y + i sin k x . On the other hand, for d+p pairing sym¬ 
metries in B i representation, we have Y 9 k = cos k x — cos k y 
and T^ = sin k y — i sin k x . The pairing amplitudes are 
related to the pairing fields by A s = 
and A p = = f( ± ^V ; cc — J^x) ( — an( i + refer to Ai and 
B\ representation respectively) with the order param- 
eters being given by if> t = \{c i+ah c^ - c i+oft c 4 ) and 

4b; = \{ci+ax tCjf + °i+axi c U )• The mean-field parame- 
ters are solved using the following self-consistent equa- 
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tions 


in = 

l^(c_ H c fet ) c °s ki, 

k 

(7) 

= 

2 ^ F^ c - fc t c fet) + (c-kic k 4 .)) sin k x , 

k 

( 8 ) 

Xoi = 

^E^) cosfe i. 

kcr 

(9) 

X11 = 

k 

( 10 ) 

5 = 

1 ~ ft y^( C lcr Cfecr )’ 

( 11 ) 


kcr 


where (c-ha c ka') and (c\ a c k a ') are computed with re¬ 
spect to the ground state of Eq.(| 6 |. Note that the 
triplet pairing can be cast into a gap function as A (k) = 
i(dk • cr)a y with dk being the d-vector. In the B\ rep¬ 
resentation, the gap function d& vector is not aligned 
with the internal magnetic field g/c. We find that the 
non-aligned d-vector lower the ground state energy in 
consistent with previous studies- 1 . When J 3 dominates, 
we find that the spin-triplet pairing wins over so that 
A 1 representation dominates. Hence as J 3 increases, su¬ 
perconducting phases will go through a transition from 
d-j-p wave (B 1 ) to s-hp wave (Hi). In the Ai representa¬ 
tion, dfc vector is parallel to vector, the ground state 
energy is further lowered down in comparison to the B 1 
representation. 



FIG. 2: Phase boundary between d + p wave (small J 3 ) to 

s+p wave (large J3). Here J10 = and J30 = — yy • Unless 
J 3 is very large, for typical strength of Rashba interaction, 
the pairing symmetry is d + p wave. 


Fig- [ 2 ] shows the phase boundary between d + p wave 
(Hi) to s + p wave (Hi) for different dopings at a fixed 
J 2 . The upper left corner of each is d + p wave while 
the lower right corner is s + p wave. It is clear that in 
the strong coupling limit, the pairing symmetry is s + p 
wave in Hi only when J 3 is at least 10 times J 30 for 
Ji = . For For typical strength of Rashba interaction, 




0.1 0.2 

lcU(eV) 


FIG. 3: Mean-field pairing amplitudes for d- wave and p- wave 
versus Rashba spin-orbit coupling strength. Here parameters 
are t = 0.45, t' = — 0.1575, t" = 0.0788 in units of eV. 


the pairing symmetry is d-\-p wave. Fig. [3] examines how 
d-wave and p-wave amplitudes change as the Rashba in¬ 
teraction aR increases. Due to that the J 3 term provides 
spin-triplet pairing channel and J 3 oc a 2 R , the p-wave am¬ 
plitude increases quadratically, A p oc a 2 R . However, for 
typical strength of the amplitude of p-wave is still 
one order of magnitude less that that of d-wave. On the 
other hand, the amplitude of d-wave also increases lin¬ 
early with aR due to the coupling of d-wave to p-wave 
via the J 2 term, which is proportional to aR. As a result, 
d-wave component is always larger than that of p-wave 
component. Hence for realistic parameter regimes, the 
resulting superconducting state is still dominated by d- 
wave component and the ground state is still in gapless 
phase. As an example, we show the boundary between 
gapped phase and gapless phase for the largest reported 
value of aR in Fig. [4j It is seen that for realistic param¬ 
eters A d > A p , the superconducting phase is gapless. 



►—►a R =0.05 

♦ ♦ a R =0.10 
*—* a R =0.15 

cx R =0.20 
■—■ a R =0.25 

• • a R =0.30 


FIG. 4: Typical boundaries between gapped and gapless su¬ 
perconducting phases with aR = 0.05 — 0.3 (a is the lattice 
constant). Here p is fixed at —0.9. It is seen that for realistic 
parameters A a > A p , the superconducting phase is gapless. 
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III. TOPOLOGICAL GAPLESS PHASES AND 
EDGE STATES 

In this section, we investigate the gapless regime of d+p 
wave. This is the regime that typical strength of Rashba 
interaction lies in. It is shown that nodal points split and 
the resulting superconducting state always supports edge 
states for any interfaces that are not exactly aligned in 
(100) or (010) directions. 


A. Bulk nodal points and their transitions in 
structure 

We start from the general bulk Hamiltonian with 
d + p mixed pairing symmetry and Rashba interaction 
in the Nambu basis H = 4^Mk)^k, where \j>k = 

(ckt? c _kf c -k4.) T - The 4x4 matrix form of h( k) is 
given by 


H k) 


( £k a R A k A p Al 'ipk ^ 

<*rA% 6fe _ 6/c ~A p A k 

A p A k -ipi c —£k o^rA p 

\ 6/c -A p Al a R A k —£ k ) 


( 12 ) 


where the lattice constant a is set to be 1, = 

sin k y + isink x , the kinetic energy is £ k = —2t(cosk x + 
cosky ) — 4t' cos k x cos k y — 2t"(cos2k x + cos2 k y ) — /i, 
and ip k is the singlet-pairing d-wave gap function. The 
gap function can be generally combined as A(fc) = 
i(ip k + d k • (j)(J y with 'ip k = —Ad(cosk x — cos k y ) and 
d k = — A p (sinfcy, sin k x , 0). The energy spectra E k can 
be solved and are given by 


Ek = ±^+^ + gl + dl±2^Al + B2, (13) 

where g fc = a R (smk y ,-smk x ), A fc = £ k g k +'ip k d k and 
Bfc = gfe x dfc. 

l«Rl>|Apl |a R |=|A p | |a R |<|Ap| 





FIG. 5: Nodal point structure in the presence of Rashba spin- 
orbit interaction. Here solid circles are nodal points for |a#| > 

|A p |, open circles are nodal points for |or| < |A p | and when 
| oh | = |A P |, nodal points in each quadrant recombine. 


In the absence of the Rashba spin-orbit interaction 
and the p-wave, the Hamiltonian do(k) supports 4 Dirac 


points determined £ k = 0 and pj k = 0. The k • p expan¬ 
sion of h( k) about one Dirac point k^> can be generally 
written as 

3 

h 0 (k) = Ywnii, ( 14 ) 

i =1 

where q = k — k^ and coordinate axes are along diago¬ 
nal axes with qi = q x + q y and q 2 = q x ~ Qy For kp in 
the first quadrant, one finds that k^> = (k Fx , k Fy , 0) = 
(kD x ,kr) x , 0). In the diagonal coordinates, we have 
v F = (2t sin k Dx + ( 2 1' + At") sm2k Dx , *-A d s\nk Dx , 0), 
7 i = s z , and 72 = s y a y . Here s is the pseu¬ 

dospin vector in the particle-hole space. In anal¬ 
ogy to the characterization of 3D Weyl semimetald 21 
and Weyl fermions in ferromagnetic superconductor^^, 
71 and 72 can be extended to be a complete set 
of 7 matrices in SU(4) representation {/, 7 ^, 7 ^} with 

{71,72,73,74,75} = {S Z ,Sy(Ty,Sy(T X ,S X ,Sy<T Z }. If WC T&- 

define operations of time-reversal and inversion sym¬ 
metries as operators with respect to q for small g, 
the Rashba spin-orbit interaction and p -wave order pa¬ 
rameter break time-reversal symmetry (operator repre¬ 
sented by i(j y K) and inversion symmetry (operator rep¬ 
resented by s x ). In this case, perturbations due to 
broken time-reversal and inversion symmetries can be 
generally expressed in terms of 7 ^ matrices defined by 
7 ij = —i/2 [j i: 7 ^pH. We find that 

h( k) = h 0 ( k) + u • b + v • b', (15) 

where b = ( 714 , 724 , 734 ), b' = ( 715 , 725 , 735 ), u = 
-smk Dx (-A p ,a R ,0), and v = sin k Dx (-A p , a R , 0). 
The corresponding energy spectrum E k of Eq. ([15]) be¬ 
comes E q and is given by 

E± = ±\Jp\ +p\ + 2 [ol\ + A 2 ) sin 2 k Dx ± 

2^2(a 2 R p 2 + A^pl) sin 2 k Dx + 4<a^A^ sin 4 k Dx , 

(16) 

where we have defined pi = v Fl qi and P2 = v F2 q2 . 
Clearly, when q\ = g 2 = 0, no longer vanishes. 
Hence the original nodal points are lifted. Similar to 
the 3D Weyl fermions, the main effect of the Rashba 
interaction and the p- wave is to split each Dirac point 
into two Weyl nodal points in 2D. Specifically, h( k) ob¬ 
tained in Eq. JT5l) corresponds to the case u || v for 3D 
Weyl semimetaPH We find that for \a R \ > |A P |, E q 

vanishes at p\ = + ^2{a 2 R — A 2 ) sin 2 k Fx and p 2 = 0; 
while for \a R \ < |A p |, E q vanishes at pi =0 and 
P 2 = ±\/2(A2-a^)sin 2 k+. 

The transformation of nodal point structure is illus¬ 
trated in Fig. [5] For a fixed a R , as A p increases, a tran¬ 
sition of nodal point configuration occurs at \a R \ = \A p \ 
before the ground state becomes a full gapped TRI su¬ 
perconducting state. The transition of nodal point struc¬ 
ture results from competition between d-wave and p-wave 
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pairing symmetries. For small p-wave amplitudes when 
|or| > |A P |, the superconducting state is dominated by 
the d-wave symmetry. Hence nodal points are determined 
by nodal lines of d -wave order parameter, which are along 
axes of k y = ±k x . The intersections of splitting Fermi 
surfaces by qlr and nodal lines of d-wave order parame¬ 
ter determines nodal points along k y = ±k x . For large 
p-wave amplitudes when |or| < |A P |, nodal points are 
determined by nodal lines of p-wave order parameter, 
which are k x = 0 or k y = 0. As a result, nodal points 
are along the qifa) axis in perpendicular to the axes of 

ky k x (ky kx') . 


B. Bulk nodes and edge states 

The nodal points revealed in the last subsection carry 
winding numbers. The associated winding numbers with 
nodes labeled by ± are shown in Fig. [6| When posi¬ 
tions of nodal points change, the distribution of wind¬ 
ing numbers also changes. It is known that positions 
of nodal points determine edge states. According to 
RefsQj and f23j for 3D bulk states, as long as the pro¬ 
jection of opposite-winding-number nodal manifolds does 
not completely overlap in the boundary Brillouin zone, 
there will be gapless surface bound states. For 2D bulk 
states, the gapless surface bound states become zero- 
energy edge states. Therefore, for a given edge, the pro¬ 
jection of nodal point to the edge determines zero-energy 
edge states. 



FIG. 6: Projection of nodal points with winding numbers la¬ 
belled by =b for (a) |or| > |A P | and (b) |or| < |A P |. Here 
each nodal point is associated with a branch cut in —k y di¬ 
rection with phase jump being ±27r. Due to the cancellation 
of phase jumps, only lines that connect ± and — are branch 
cuts with non-vanishing phase jumps. The projection of each 
line connecting ± and — onto the edge V axis gives rise to one 
zero-energy flat band at the projected range of wavevector p y . 
Absolute values of numbers (±2 or ±1) shown in the edge T 
axis are numbers of flat-bands in the corresponding p y range. 
It is seen that for realistic Rashba strength, |qr| > |A P |, one 
always gets single zero-energy mode for certain wavevector p y 
along the edge T. 


Fig- © illustrates the projection of nodal points to 
edges. For each nodal point, one associates a branch cut 
with phase jump being ±27t. In Fig. [6j branch cuts are 
chosen in —k y direction. Due to the cancellation of phase 
jumps, only lines that connect ± and — are branch cuts 
with non-vanishing phase jumps. Here the phase jumps 
for the line connecting ± to — and that for — to ± are 
opposite in sign. By projecting branch cuts onto the mo¬ 
mentum line (denoted by p y ) that represents the edge, it 
is clear that for \cxr\ > \A P \ shown in Fig. [6h, the nodal 
manifolds for ± nodes do not overlap completely except 
when orientations of edges are exactly along (100) or 
(010) directions. Therefore, for p y in the non-over lapped 
region, there always exists single zero-energy mode along 
the edge. The single zero-energy mode is a Majorana 
mode. More specifically, the projection of each connec¬ 
tion of ± and — to a given edge Y gives rise to one zero- 
energy mode for each wavevector p y along Y that lies in 
the non-over lapped region. It thus results in a flat band 
at zero-energy^. Therefore, in the overlapped region, 
there will be two flat-bands; while in the non-over lapped 
region, there is only one branch flat-band which is a Ma¬ 
jorana flat-band. As a result, we find that except when 
the edge is exactly along (100) or (010) directions, in gen¬ 
eral, there always exist Majorana modes for any edges. 
Since real interfaces always possess facets in different ori¬ 
entations, existence of Majorana modes appear to be a 
robust feature. 

As the amplitude of p-wave increases until |A P | > |or|, 
nodal points switch to be in the k 2 direction. As illus¬ 
trated in Fig. [6)3, the nodal manifolds for ± nodes do not 
overlap completely except when orientations of edges are 
exactly along (110) or (110) directions. Therefore, except 
when the edge is exactly along (110) or (110) directions, 
in general, there always exists Majorana modes for any 
edges. Note that in the case when |A P | > |or|, there 
is a region near p y = 0 for (110) edge in which as we 
shall show in below, the connection of nodal points with 
the same winding number leads to dispersive Majorana 
modes. 


C. Edge states along (110) edge 

In this subsection, we shall verify connections of edge 
states to the nodal structure for the particular orientation 
of edge (110). The (110) edge is known to be the most 
important edge that exhibits zero-bias peak in the tun¬ 
neling spectrum due to the existence of Andreev bound 
states at zero energ^PH. It is therefore important to ex¬ 
amine how edge states change when the Rashba spin- 
orbital interaction is included. We shall start by partial 
Fourier transforming the Hamiltonian along the (110) in¬ 
terface direction. Let p y be the Fourier wavevector along 
(110) direction. The Hamiltonian can be expressed as 
H = J2 Py (Py) Hid (Py)^(Py) and is characterized by 
an one-dimensional Hamiltonian HmiPy) for a given p y , 
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which is given by 

Hi d(p v ) 


( h(p y ) Mp V ) \ 

V M Py ) ~ hT (-Py) j ’ 


(17) 


where the hopping and pairing matrices h(p y ) and A (p y ) 
are given by 


k(Py) = 


htt hU 

ri Py ri Py 

htf hH 

Py Py 


. MPy) = 


A tt A n 

Py Py 


(18) 


Py 


Py 


The block matrices h(p y ) are expressed as 


( t_ l ijj ij| 0 ... ^ 



*11 * L *11 *|| '• 

*|| *11 * J - *11 *|| 

0 if, t|| i± t„ 



\ ; o t' t|| / 


( 0 

t\ 

0 

...\ 

ti 

0 

ti 

0 

0 

h 

0 

h 

V : 

0 

t2 

‘•J 



with and /i** = (/i);() t. The pairing matrices 


Py 


Py 


Py ' 


(d) r|=0.5 


Py 




0 

A 2 

0 


A! 

0 

a 2 




0 

Ai 0 
0 Ai 


A n = 

’ Py 


: 0 A2 '•• / 



(e) r|=0.36 



Here elements in matrices are given by ij_ = 
-2t'cos(p y S' 2y ) - /i, £|| = -2tcos(p 2/ ^i 2/ ), tj| = -t' - 

2t" cos(p y 5i y ), ti = ^(-ie^y^y + e-^ 1 **), t 2 = 
oLKfa-WySiy _ e % Py^iy), Ai = — A 2 = iA^sin^d)^), 
A 3 = ^(— ie ip y 6l y — e^Py^y), A 4 = ^^ie~ ip y 6l y + 
e l P y Siy ) with = 1/2 and 5^ = ^ = 1- 

We shall focus on the case when the superconducting 
side is semi-infinite. For a given p y , the system becomes 
one-dimensional and its topology can be characterized 
by the winding numbed 25 . Here we extend the winding 
number calculatiorP^ to any edge that makes an angle 6 
clock-wisely with respect to (100) direction. The winding 
number W is then given by 

1 r 

W(p\\,0) = —Tra J dp ± d p± lndetg(p||,pj_,0), 

1 N 

= 2 tt / d P ±eabTn *(P\\’P±’ e } d P-L m b(P\\’P-i’ e )- 

(19) 

Here p\\ is the wavevector along the edge, p± is the 
wavevector perpendicular to the edge, and the matrix q is 
given by q = *(£(p, 0)-iip( p, 0))a y + i(g(p, 9)-id( p, 0)) ■ 
aery. mi(p,6 > ) and m 2 (p,0) are the real and imaginary 
parts of the phase factor of det q( p, 6) respectively. 

For a give p y , if number of sites in perpendicular to 
(110) edge is TV, H\o{p y ) is a 4TV x 4 N matrix and can 


FIG. 7: (a) Comparison and consistent check of analysis by 
projection and analysis by winding number for ( 110 ) edge. 
(b)-(e) Evolution of quasi-particle spectrum and edge states 
for ( 110 ) edge of cuprates in the gapless phase when aR 
changes. Here 77 = olr/ A p . The mean-field solution corre¬ 
sponds to 77 = 12.86 with S = 0.16 in (b). Edge states are 
labelled as blue and red solid lines with red lines representing 
Majorana edge states. It is clear that Majorana flat-bands at 
zero energy exist only for physical accessible region aR > A p ; 
while the dispersive Majorana modes exist near p y — 0 only 
for aR < A p . 


be exactly diagonalized. For large TV, one obtains numer¬ 
ical solutions for semi-infinite superconducting state. In 
Fig. [7Fa), we show how the quasi-particle spectrum for 
the (110) edge changes as the Rashba spin-orbit inter¬ 
action aR varies. Starting from the mean-field solution 
with 77 = aR/A p = 12.86 shown in Fig. [T^a), there are 
two Majorana flat-bands with winding number \W\ = 1 
at 2 < p y < 2.6. The decreasing of aR leads to the 
reduction of the Majorana flat-band regimes as shown 
Fig.JT^b). In particular, as indicated in Fig.^c) and (d), 
Majorana flat band disappears for 77 < 1 (i.e. aR < A p ). 
Instead, we see that dispersive edge states occur near 
p y = 0 with 77 < 1 . 

To reveal the nature of the dispersive edge states, we 
fix aR and varies A p . As illustrated in Fig. [7] (a)-(c), 
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FIG. 8 : Illustration of the nature of dispersive edge states for 
(110) edge in the gapless phase. Here t = 1, t' = t" — 0, /i = 
0.5, and A v > aR with A a *= 0.1 and otR = 0.05 being holding 
fixed. As shown in (a)-(c), as Ap increases (A p = 0.08, 0.2 
and 0.28), =b nodal points from different quadrants move to¬ 
ward each other and annihilate. Correspondingly, as shown in 
(d)-(f), the region with dispersive edge modes around p y = 0 
expands and the region (labelled by d=2) with fermionic An¬ 
dreev bound states shrinks. It is clear that dispersive edge 
modes near p y = 0 evolve into helical Majorana modes of the 
fully gapped p -wave superconductors. 


when the amplitude A p increases, db nodal points from 
different quadrants move toward each other and annihi¬ 
late. Correspondingly, as shown in (d)-(f), the region 
with dispersive edge modes around p y = 0 expands and 
the region with fermionic Andreev bound states shrinks. 
It is clear that dispersive edge modes near p y = 0 evolve 
into helical Majorana modes of the fully gapped p-wave 
superconductors. Therefore, when A p > a^, dispersive 
edge states around p y = 0 emerge and one expects that 
they are of the same nature as helical Majorana modes 
of the fully gapped p-wave superconductors. 

To examine if these dispersive edge mode are in¬ 
deed Majorana modes, one needs to examine if the 
quasi-particle obeys rf(p y ,E) = 7 (~p y , -Ef- 1. In gen¬ 
eral, the eigenstate at i site to the one-dimensional 
Hamiltonian Hm(p y ) is represented by a 4-component 
(^1 (p y , E, i) , v 2 (p y , E, 1 ) , Ml (py, E,i) , u 2 (p y , E, i) ) so that 


the quasi-particle is given by 

7 \p y , E, i) = Uxipy, E ,+ U 2 (p y , E, i)cj yi4 

+v 1 (p y , E, i)c- Py rf + v 2 {p y , E, i)c- Pytii . (20) 

For Majorana modes, one requires 7 ^(p y ,E,i) = 
7 {~Pyi—E,i), which in turn requires Vi(p y ,E,i) = 
u i and v 2 (p y ,E,i) = u$(-p y , -E,i). How¬ 
ever, due to an uncertainty in defining Ui and Vi up to a 
U( 1) global phase factor, in order for 7 to be Majorana 
modes, one requires 


Vi (p y ,E,i) = e l ^u\{—p y , —E,i), 
v 2 (p y ,E,i) = e'*i4(-py,-E,i). (21) 

Therefore, \vi/ui \ = \v 2 /u 2 \ = 1 is the condition for edge 
modes to be Majorana modes. In Fig.[9j we examine en¬ 
ergy dispersions for (110) edge near p y = 0. As shown in 
the inset, \vi/ui \ = \v 2 /u 2 \ = 1 is satisfied for dispersive 
edge modes. Hence when A p > the emerging disper¬ 
sive edge states around p y = 0 are dispersive Majorana 
modes. 



§0.8- 

— 0.6 ^. 

0 0.05 0.1 0.15 0.2 0.25 

b 

FIG. 9: Energy dispersions for (110) edge near p y = 0. Here 
solid lines represent bulk states. Filled circles represent edge 
states. Lower panel: The Majorana condition, 7 ^(E,p y ) = 
7 (—E,—p y ), is checked by examining \vi/ui\ = \v 2 /u 2 \ = 1. 
Clearly, the dispersive edge states are Majorana modes. 


IV. TUNNELING SPECTROSCOPY 

The tunneling spectroscopy is one of means to examine 
Majorana fermions. A typical experiment is the measure¬ 
ment of the tunneling spectrum for normal metal (V), in¬ 
sulator (/) and superconductor (S) junction, i.e., the NIS 
junction, in which zero-energy edge states would appear 
as a zero-bias conductance peak. However, the zero-bias 
conductance peak may also arise from other mechanism 6 
and hence it is not considered as a smoking-gun evidence 
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for Majorana fermions. On the other hand, the 47r peri¬ 
odicity of super current flowing across the Josephson junc¬ 
tion is considered as a smoking-gun evidence for existence 
of Majorana fermions. In this section, we shall examine 
both the spectrum of the NIS junction and Josephson 
junction for high T c cuprates in the presence of Rashba 
spin-orbit interaction. It will be shown that dispersive 
Majorana modes would exhibit a unique feature as a 
small plateau in tunneling spectrum near zero bias peak. 
Furthermore, we find that the flat-band Majorana modes 
always result in 4tt -Josephson effect in typical Joseph¬ 
son junctions. Therefore, by suitable designing junctions 
involved in the tricrystal experiments, the flux trapped 
in the center of the tricrystal jumps in unit of two flux 
quanta, which would be a convincing way to detect Ma¬ 
jorana fermions. 


by fermionic and Majorana edge states. Both of them 
contribute to the zero-bias peak based in the Blonder- 
Tinkham-Klapwijk theory 30 . On the other hand, the 
dispersive Majorana fermions occur at finite energy scale 
and hence its appearance would be unambiguous. 

For a planar NIS junction, after partial Fourier trans¬ 
formation along the edge of the junction with p y being the 
wavevector along the edge, the junction is effectively a 
one-dimensional system. For a give p y , the effective one¬ 
dimensional Hamiltonian for the superconducting side 
with (110) edge is given by Eq.(17), HioiPy)- Follow¬ 
ing Blonder, Tinkham, and Klapwijk 30 -!, the NIS junction 
is formed by connecting Hid{p v ) with the effective ID 
Hamiltonian of the normal metal side. If the one di¬ 
mensional eigenstate T(p y ) at j site is represented by a 
4-component ifjj . The tight-binding Bogoliubov equation 
can be generally written as 


A. Tunneling spectrum of NIS junction 


The well-known zero-bias conductance peak observed 
in high-T c cuprates along (110) direction manifests the 
fermionic^ Andreev bound states at the i nterf ace of NIS 
junctional As indicated in the Sec. Ill C, in the 


presence of the Rashba spin-orbit interaction in high- 
Tc cuprates, (110) edge also hosts Majorana flat bands 
when p-wave is small. The detection of this Majorana 
flat bands is thus ambiguous in the NIS planar junctions 
due to the coexistence of zero-energy flat bands formed 


T^j-i + Tt^+i + T'^j- 2 + T'^j^ 2 + fi'ipj 

+V6 Xji - 1 if j = E4> j , (22) 

where the potential V = Vr z <S> (Jo is the tunneling bar¬ 
rier at the interface due to the insulator layer, the su¬ 
perconductor S occupies j > 0 and the normal metal N 
occupies j < 0. Matrices T and T' denotes effective near¬ 
est neighbor and next nearest neighbor hopping, while fa 
represents the effective on-site potential. They are given 
by 



h 


A t i 

A s 

\ 


( 

t± 

0 

0 

0 

\ 


(tl 

0 

0 

0 N 



h 

-A, 

A(2 


,fa = 


0 

t± 

0 

0 


,T' = 

0 

'll 

0 

0 


— A 4 2 

—A s -t|| 

t\ 



0 

0 

-t± 

0 


0 

0 

-'ll 

0 


A s 

—Ati 

t2 


J 


V 

0 

0 

0 

-t _L 

J 

1° 

0 

0 

-tl / 


Here for the metal side, corresponding matrices are de¬ 
noted by Tn and fax • A s , Ati, A*2, and £j| are set to 
zero, while ty, t_ l, £i, and t 2 are set to be hopping ampli¬ 
tudes -2t N cos(p y 5 ly ), —/in , %f(- ie ip y Sl y + e~ ip y 6l y), 
and ^[ie~ ip y 8l y — e iPySly ). For the superconducting 
side, corresponding matrices are denoted by Tg, fas and 
T' s with £|| = —2tco$(p y Siy), tj| = —t' — 2t" cos(p y 5i y ), 
t _L = -2t'cos(p y S' 2y ) -//,*! = °f(-ie ip y 5 'y + e ~ ip y^y), 
t 2 = ~f(ie~ ip y 5l y - e ip y 5l y ), A s = -iA d sm(p y S ly ), 
A n = % ’-(ie-iPyhy+jPyhy ) and A t2 = %(-ie“^ <5l ^ + 

giPySly ^ . 

At a given energy E, far away from the interface, the 
quasi-particle wavefunction ifjj satisfies the bulk Bogoli¬ 
ubov equation. As illustrated in Fig. [lOj there are two 
incident waves (indicated as black and blue dots) in the N 
side, which will be denoted by 4^ and 4>g 2 . For reflected 


waves, there are 2 bulk quasi-particle wavefunctions indi¬ 
cated by black and blue dots and 2 bulk quasi-hole wave- 
functions indicated by open circles. The corresponding 
normalized wavefunctions are denoted by 4> ei , 4> e2 , 4^, 
and 4>h 2 . Similarly, for the S side, there are 2 particle-like 
and 2 hole-like wavefunctions denoted as 0 ei , 0 e2 , <^ 15 
and (j)h 2 - By using these wavefunctions, i/j j(p y ) can be 
expressed as 

for j > 0, 1pj{Py) = tnZi4>n(Py,Zn), 

n=ei,e2,hi,Ji2 

for j < 0 , i>j(p v ) = ^j(p y ) 

= ^ ' T n Zn*^n(Pyi Z n ). (23) 

n=ei,e2,ei,e2,hi,h2 

Here indices with e\ and e 2 represent quantities of in¬ 
cident waves. The wavevectors p™ and p™ in the corre- 
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sponding quantities z n = e w ^ and z n = e ip * are per¬ 
pendicular to the interface in the 5 side and the N side 
respectively, p x and p x must be solved from the bulk en¬ 
ergy dispersion E(p™,p y ) = E. t n and r n are correspond¬ 
ing transmission and reflection amplitudes for c/) n (p y , z n ) 
and & n (p y , z n ). Solutions for j > 0 and j < 0 with forms 
given by Eq.(23) are then joined together at the interface 
j = 0. By comparing Eq.(22) at j 0 and j <C 0 with 

Eq.P 


at j = — 1 and j = 0, we find that the boundary 
conditions are given by 


FIG. 10: Schematic diagrams of the NIS junction for quasi¬ 
particles traveling (right arrows) from N to S and reflected 
(left arrows) from S to N. Here upper panel is for otR > A p 
and lower panel is for aR < A p . In the N side, solid lines and 
dash lines represent spectra of particles and holes respectively. 
Different colors (black and blue) stand for two bands split 
by the Rashba interaction. Black and blue dots with arrows 
represent incident and reflected particles, while open circles 
denote the reflected holes. In the S side for E > Eq (p y ), 
gray dots stand for the particle-like and open circles denote 
the hole-like transmitted modes. For aR < A p , only two 
nodal points are present in S side. Specifically, two evanescent 
modes represented by dash arrows can be found by solving 
E(Px,Py) — E with E > Eo(p y ). Here, E(p™,p y ) is the lowest 
positive E p in Eq.(13). 


- cdr 1 ^-!, 

(24) 


After re-arrangements of the above two boundary condi¬ 
tions, amplitudes t n and r n can be solved by the following 
8x8 linear equation 


^S Z e 2 1( t ) e2 


[J e 2 

B ~-i 


0/ii 


<K 2 


U ei <S> ei 

-Tlz~^ ei 


^e2 thi ^h 2 p e\ p e 2 p h\ p h 2 ) 


Ue 2 ®e 2 

~^'N Z e 2 ^e 2 

-Ue i^ei - 




- lA e2 *&e 2 


Uh 2 $h 2 


^N Z h 2 ^ h 2 


i +Z E ^®e 2) 


(25) 


where U n = (T^) 1 Vz n 1 - 1. 

For a given p y , there is a minimum of bulk energy 
at Eo(p y ) as illustrated in Fig. [lO] When the energy 
E of an injected particle with E > Eo(p y ), p x and p x 
are both real. We obtain two transmitted electron-like 
and two hole-like quasiparticles in the S side and two 
reflected electrons and holes in the N side. On the other 
hand, for E < Eo(p y ), two electron-like and two hole-like 
in the S side become evanescent modes with p e x and p x 
being complex numbers. Note that for aR < A p , even for 
E > Eo(p y ), since there are only two nodal points, both 
transmitted and evanescent modes are present in the S 
side. 


Solving Eq. (25), one obtains solutions for t n and r n 


which then determine the transmitted and reflected par¬ 
ticle currents given by 


Jj,j' ^ 


i (BBb - tyM?) , 


(26) 


in the N side and similar expression for the S side with 


T/v and \F/ replaced by Ts and ijjj. Hence we obtain the 
particle current J t e through the tunneling barrier V in 
the N side and the reflected particle and hole currents: 

and J x . The normalized tunneling conductance for 
the NIS junction is then given by 30 

vn(E,V) = * y n (E,V,p v ), (27) 

Py 

where g 0 (E, V) = J2 Py Jt( E , V ’Py)/ J i( E >Py) is the bal- 
listic conductance of the N side and <j n (E,V,p y ) is de- 
fined to be 1 - [J?(E,V,p y ) - Jp(E,V,p y )]/J?(E,p y ). 

In Fig. 13 we show the computed tunneling conduc¬ 
tance. In the strong interface scattering limit V —> oc, 
the tunneling conductance isproportional to local density 
of state po(E) at (110) edge 3 -^, dominated by edge states. 
For edge states with the dispersion E(jp y ) oc p™, the local 

density state is po(E) oc E ~^ as illustrated in the in¬ 
set of Fig. 13 Clearly, the zero-bias conductance peak in 
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Fig. [TTImanifests the dispersionless Andreev bound states 
shown in the inset of Fig. [IT] while the small plateau at 
finite energy scale reveals the linearly dispersive Majo- 
rana fermion near p y = 0 . As a consequence, an overall 
candleholder-like structure of the tunneling conductance 
reflects the coexistence of Andreev flat-band and disper¬ 
sive Majorana fermions at the interface. 



FIG. 11: Tunneling conductance near zero bias for Metal- 
Insulator-high T c cuprate superconductor in the presence of 
Rashba spin-orbit interaction (qlr < A p ). Here E is measured 
in terms of t with t' — t" — 0 ,aR — 0.13, A^ — 0.32, A p — 
0.16. The tunneling conductance <jn near E — 0 exhibits 
a candleholder-like structure, which reflects the local density 
of states at the edge in the strong interface scattering limit 
(V —oo). The insets show the corresponding local density 
of states and the energy spectrum in which the Andreev flat 
edge modes contributes to the zero-energy divergence while 
the dispersive edge modes contribute to the plateau. 


B. SIS 7 -Josephson Junction 

In this subsection, based on our mean-field solutions, 
the quasi-particle spectrum and Josephson current in 
SIS 7 junctions are investigated. We shall examine pe¬ 
riodicity of Josephson current in SIS' junctions with S 
and S' formed by high T c cuprates in different edge ori¬ 
entations as illustrated in Fig. [l3|^a). In particular, we 
will show that the hybridization of Majorana zero-energy 
flat-bands of S and S' leads to 47 r-periodicitjW. How¬ 
ever, we show that this only occurs when there is an 
overlap between Majorana zero-energy flat-bands from S 
and S'. Based on these findings, we explore and design 
the tricrystal junction which can generate a spontaneous 
vortex with half-flux quanta with jump in unit of two flux 
quanta as the magnetic field is applied. This could serve 
as a way for detecting the elusive Majorana fermions in 
the high-T c cuprate-based heterostructure. 

We start by representing the wavevector along the 
junction by p y . For a given p y , the thin insulating layer 


/ and superconducting sides S and S' are effectively one¬ 
dimensional with coordinates denoted by x and their cor¬ 
responding Hamiltonians are generally given by Eq. 
denoted by Hi(jp y ), Efs(p y ), and Hs f (p y ) respectively. 
The effective Hamiltonian of SIS' is then given by 

HsiS'iPy) = H s {p v ) + H^Py) + H S '{Py) 

+ Hsi{p y ) + His'(p y ) + H.C., (28) 

where the insulator occupies — xjy < x < xn , S oc¬ 
cupies x < — xat, S' occupies x > x^, and Hsi(p y ) 
and His'(p y ) are the tunneling Hamiltonians describing 
tunneling between superconductors and the thin insula¬ 
tor. To investigate the Josephson current, we impose a 
a U(l) phase difference 0 = cj)s' — 4>s = 4> between the 
pairing potentials in S and S' sides by setting the pair¬ 
ing potential = A (p y )0(—Xi — xn) for S side and 
A s f = A (p y )e l ^0(xi — xn) for S' side. 

(a) (b) 

Py 



FIG. 12: (a) Schematic plot of a SIS 7 junction. Here the rela¬ 
tive angles of (100) directions (denoted by a and a' axes) of S 
and S 7 to the tunneling direction (x axis) are 0 and 0' respec¬ 
tively. (b) Comparison of projected p y along the junction for 
edge states that would appear in isolated and semi-infinite S 
and S 7 . Here 6' > 0 and red and blue dots are nodal points. It 
is seen that p y s of Majorana zero-energy modes (flat-bands) 
from S and S 7 may overlap in the region p 2 < \p y \ < pi and 
result in 47r periodicity in Josephson current. 


We first analyze edge states of isolated semi-infinite su¬ 
perconducting sides (S and S 7 ). Following Fig. [ 6 } by pro¬ 
jecting nodal structures onto an edge, one obtains possi¬ 
ble number of edge modes for each p y along the edge. In 
Fig.^b), we show the comparison of two edges with dif¬ 
ference orientations. It is seen that a junction formed by 
these two edges may hybridize two Majorana zero-energy 
flat-bands (e.g. P 2 < \p y \ < pi), one fermionic Andreev 
flat-band with one Majorana zero-energy flat-band (e.g. 
Pa < \Py\ < P 3 ), or two fermionic Andreev flat-bands 
(e.g. \p y \ < £> 4 ). As we shall see, only the overlap of zero- 
energy Majorana bound states lead to their hybridization 
and result in 47 r periodicity of Josephson current when 
these two edges form a SIS 7 junction 34 . If the relative 
angle of ( 100 ) directions of a superconductor S to the 
tunneling direction (x axis) is #, the relative angle of pro¬ 
jected p y to the (100) direction is 7 r /2 — 0. Hence projec¬ 
tion of the Majorana zero-energy flat-band onto the junc¬ 
tion interface is k~ cos( 7 r /4 — 0) < p y < cos( 7 r /4 — 0), 
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where k± are magnitudes of wavevectors for nodal points. 
Assuming that both S and S' have the same doping level, 
the projection of Majorana zero-energy flat-band for S' is 
also determined by k ± and hence the projected Majorana 
zero-energy flat-band is given by k~ cos(7r/4 — O') < p y < 
k + cos(7r/4— O'). For doping S ~ 0.16 and qlr ~ 0.05—0.3, 
we find that k~/k + ~ 0.92 — 0.77 and the Majorana zero- 
energy flat-bands for the SIS' junction overlap as long as 

\0 — 0'\ <21° -39°. (29) 


To illustrate effects of overlapped zero-energy Majo¬ 
rana flat-bands in the junction, we examine quasi-particle 
spectrum and Josephson current for junctions : (i) 0 = 
45°/(110) edge, 0' = 0°/(100) e dge and (ii) 0 = 45°, 


13'a), we show quasi- 


0' = tan -1 2/(210) edge. In Fig. 

particle spectrum for the S(45°)-I^S r (0°) junction. Here 
due to finite lengths of S ans S' (Nx = 625), in-gap 
states include edge states from outer edges of S and S'. 
It is seen that there is no Majorana zero-energy state 
(single branch) formed by hybridization of edge states 
from S ans S'. This is verified by checking energies of 
edge states versus the U( 1) phase difference <p shown in 
Fig. [l3tb). The period is 27r. Next, we consider the 


(a) 



p =-2.5035 


0.06 

0.04 

0.02 

! ° 
-0.02 
-0.04 
-0.06 


<t>fa) 


(b) 




0.04 

0.02 

0- 


-0.02 

-0.04 



FIG. 13: (a) Low-lying quasi-particle spectrum for S(45°)-I- 

S'(0°) junction with <p = 0.57T. Here due to finite lengths of 
S ans S' (Nx = 500), in-gap states include edge states from 
outer edges of S and S'. Edge states inside the junction are 
marked by blue (Andreev Fermionic modes) and red (Majo¬ 
rana modes) lines. It is seen that single branch Majorana 
flat-band is missing, (b) Energies versus the phase difference 
<p for edge states in the junction for a few p y ’s. The periods 
are found to be 2n. Here length of S and S' is N x — 625. 


S(45°)-I-S'(tan 1 2) junction. Fig. 14 'a) shows the full 


quasi-particle spectrum. It is seen that there is a Majo¬ 
rana flat-band for 2 < \p y \ <3, resulting from hybridiz¬ 
ing two Majorana zero-energy flat-bands from S and S'. 
In addition, there are in-gap states due to edge states 
in the junction formed by hybridizing one fermionic An¬ 
dreev flat-band with one Majorana zero-energy flat-band 
or two fermionic Andreev flat-bands. In Fig. 14 T), we 


examine energy versus phase difference <p for edge states 
in the junction. It is clearly seen that only Majorana 
modes due to hybridization of Majorana fermions from S 
and S' exhibit 47r periodicity. The 47 t period results from 
conservation of the fermion parity^ 7 and the fermion 
parity is conserved only when two Majorana zero-energy 


modes are hybridized. We thus verified that only when 
the Majorana zero-energy flat-bands for each edge that 
form the junction overlap, the supercurrent in the junc¬ 
tion shows 47r periodicity. 



FIG. 14: (a) Quasi-particle spectrum of S(45°)-I-S'(tan _1 2) 
junction for <p = 0.57T. Here length of S and S' is N x = 500. 
Edge states inside the junction are marked by blue (Andreev 
fermionic modes) and red (Majorana modes) lines, (b) Energy 
versus phase difference <p for edge states in the junction for a 
few p y ’ s. Here length of S and S' is N x = 625. Plot of p y m 
—2.4504 represents the spectrum due to the hybridization of 
Majorana fermions from S and S '. The red and black solid 
lines denote different fermion parity P = ±1. Plots of p y = 
— 1.8535 and p y = —1.665 stand for the spectrum formed by 
hybridizing fermionic Andreev zero-energy bound states. The 
remaining plot for p y = —1.2881 is the spectrum resulted from 
hybridization of fermionic and Majorana zero energy modes. 
It is seen that only hybridization of Majorana fermions results 
in 47 t periodicity. 


The 47r periodicity implies fluxes of vortices trapped 
in the ring formed by the SIS' junction jump in unit 
of two flux quantsP. One of configurations that allows 
one to explore the pairing symmetry and also vortices 
trapped in high T c cuprates is the tricrystal experiment 
with a 7r-ring 31 . Here we examine the periodicity of 
1/(1) phase winding phase for the tricrystal configura¬ 
tion shown in Fig. [l5|a). As illustrated in Fig.p^a), for 
general tricrystal configurations, one may fix the edge di¬ 
rection of the cuprate in the bottom labelled by C to be 
in (010) direction. Let angles between a axes of the upper 
left crystal (A) and right crystal (B) and (100) direction 
of the bottom crystal be 6a and Ob • If the interface be¬ 
tween A and B be T and the angle between T and the 
(100) direction be #, the interface T is an S(0 — 6a)- I- 
S '{0 — Ob) junction. Clearly, since there is no Majorana 
zero-energy mode in the (010) edge, in the tricrystal con¬ 
figuration shown in Fig. [f~5^a), Majorana modes can ap¬ 
pear only in the interface between A and B. According 
to Eq.(29), as long as 0 ^ Ob and 0 ^ 6a, the condition 


\\0a-0\-\0 b -0\\ <21°-39° 


(30) 


is satisfied, there are Majorana zero-energy modes in the 
interface between A and B. On the other hand, the su¬ 
percurrent passing an junction interface between super¬ 
conductors i and j is given by I = I c cos 20i cos 20j (clean 
limit) or I = / c eos2(6^ + 0/) (dirty limitin which 
Oi and Oj are angles of the crystallographic axes (such as 
(100)) with respect to the junction interface. For the ring 
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FIG. 15: (a) Configuration of tricrystal experiment. Here 6 a 
and 6b are the angles between the a axis of the upper left 
crystal (A) and right crystal (B) and (100) direction of the 
bottom crystal, (b) Regions of parameters for 7r rings with 
Majorana modes (shaded areas with dots). Here 6 ~ 6b- 
Grey areas are regions for zero rin^ 32 l White areas are regions 
for 7r rings without Majorana modes, and shaded areas with 
lines are regions for zero rings with Majorana modes. For 
doping S = 0.16 and <ar = 0.05 — 0.3, the width W for the 
region with Majorana fermions is in the range of 42° — 78°. 


(a) (b) 



FIG. 16: (a) Energy gap A t due to the coupling of S and 
S 7 in the SIS' junction. Here 6 is fixed to be 30° and the 
energy for the quasi-particle of S is chosen at the nodal point 
with the Fourier wavevector p y . At is the minimum energy 
of quasi-particles in S 7 with fixed p y . At generally does not 
vanish except when S and S' are the same or 6 + 6' = 90°. 
(b) Comparison of the enegy-phase relation due to quasi¬ 
particles (solid circle) and Majorana modes (open circle) lo¬ 
calized at the SIS 7 junction. The Josephson current due to 
quasi-particles is heavily suppressed. 


which, when combined with with Eq.(32), determines pa¬ 
rameters that allow 7r rings to host Majorana fermions. 
In Fig.^b), we summarize regions of parameters for n 
rings and zero rings with Majorana fermions. It is seen 
that there are large parameter regions satisfying the re¬ 
quirement for hosting Majorana fermions. Therefore, one 
expects that the half-flux vortex trapped in a typical i r- 
ring would jump in unit of two flux quanta in external 
magnetic fields. 


C. Influence of quasi-particles and k z dispersion 


configuration shown in Fig. [T5^a), in order that there is 
effectively a 7r-junction in the three junctions involved, 
we require cos 2 2 6c cos 2 6a cos 20b cos 2 (6 — 6a) cos 2(6 — 
6b) < 0 (clean limit) or 008 2 ( 6 ^ + Oa) cos2(6c + 
6b ) cos 2(26 — 6a~Ob) <0 (dirty limit). Hence to form a 
half-flux vortex in the ring, in the clean limit, one requires 
cos 20 a cos 26b cos 2 (6 — 6a) cos 2 (6 — 6b ) < 0 . Therefore, 
it is sufficient to require one of angles, 6At Ob, 6 — 6a, 
and 6 — 6b to be greater than 7 t/ 4 and all the others to 
be less than 7 t/ 4. There are two possibilities: 

either 7 t/4 < 0a/b < tt/2, 

or 7 t /4 < e - e A/B < tt/2, (31) 

and all the other angles are less than 7 t/ 4. As a result, 
the 7r-junction is either the junction between A (or B) 
and C or the junction between A and B. For typical con¬ 
figurations of tricrystal experiments^, one has 6 ~ 6b- 
The requirement of a ring to be a 7r-ring for the clean 
and dirty limits becomes the same and is given by 


cos 26a cos 26b cos 2(6 a — 6b) < 0. 


(32) 


On the other hand, Eq.(30) implies \6a~6b\ < 21° — 39°, 


In this subsection, we examine stability of 47r period¬ 
icity and effects due to energy dispersion along c axis. 
Specifically, we examine whether the 47 t periodicity is 
spoiled by the presence of quasi-particles. 

We first note that the 47r periodicity relies on the con¬ 
servation of Fermion pairty. In the presence of quasi¬ 
particles, the fermion parity is no longer conserved 7 . 
However, if the energies of quasi-particles are gapped, 
the fermion parity is approximately conserved at low 
temperature^. As indicated in Fig. [IJ the bulk super¬ 
conducting states in realistic parameter regime are gap¬ 
less. However, the SIS 7 junction involves the coupling 
of the superconducting states in two different orienta¬ 
tions. When two superconducting states couple, since 
the Fourier wavevector p y that is parallel to the junc¬ 
tion is conserved, only quasi-particles with the same p y 
couple. Due to different orientations of S and S 7 relative 
to the junction, p y of nodal points for S and S 7 are not 
the same. As a result, quasi-particles at nodal points 
of S (or S 7 ) couple to quasi-particles at finite energies 
of S 7 (or S). Hence the resulting quasi-particles residing 
near the junction are gapful and there is an energy gap, 
A t, associated with each orientation of the tunneling 
junction. In Fig{l6^i, we examine A t versus 6' by fixing 
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0 = 30°. Here the energy of the quasi-particle for S side 
is chosen at the nodal point with the Fourier wavevec- 
tor p y and A t is defined to be the minimum energy of 
quasi-particles in S' with fixed p y . It is seen that ex¬ 
cept for special orientations when S and S' are the same 
or 0 + O' = 90° , At generally does not vanish. The 
parity conservation thus holds for quite general orienta¬ 
tions of S and S'. In addition to the protection of the 
pairty conservation, the momentum mismatch also sup¬ 
presses the quasi-particle contribution. This is illustrated 
in Fig. 16 3. It is seen that E(^>) for the contribution of 
quasi-particles is much more flat in comparison to that 
due to Majorana modes. Hence the Josephson current 
( dE (i6 * ) due quasi-particles is heavily suppressed and 


leads to little effects. 

Finally, we discuss effects of the energy dispersion 
(k z ) along c-axis. It is generally accepted that high- 
Tc cuprates are approximately two dimensional with 
weak dispersion along c-axis. Therefore, results de¬ 
rived in the above are approximately correct. For 
high-T c cuprates with one Cu02 per unit cell, the en¬ 
ergy dispersion along c-axis can be approximated h jm 
A£k = —2t z cos(k z /2)(cosk x — cos k y ) 2 cos k x /2 cos k y /2 
with t z ~ 0.1 1. Following the discussion on Eq. @ , for 
a given k z , nodal points are determined by £& + = 0. 

In other words, A£*. changes the Fermi surface. For real¬ 
istic parameters, \(Xr\ > A p , nodal points are along axes 
k y = ±k x on which A£/~ = 0 is satisfied. Hence nodal 
points have no k z dispersion. Therefore, all results based 
on projection of nodal points (including 47r periodicity) 
are the same except that there are more Majorana edge 
modes due to the k z dispersion. On the other hand, for 
A p > |or|, the k z dispersion changes the Fermi energy 
by 10% for a give k z . As a result, there is no qualitative 
change on the dispersive Majorana modes except for the 
corresponding 1% change of the range for the existence 
of the dispersive Majorana modes. 


V. DISCUSSION AND CONCLUSION 

In conclusion, we have explored effects of the Rashba 
spin-orbit interaction on the tunneling spectroscopy of 
high-T c cuprate superconductors. The mean-field phase 
diagram in the large-U limit of Hubbard model is es¬ 
tablished. It is shown that due to the Dzyaloshinskii- 
Moriya and spin dipole-dipole interactions induced by the 
spin-orbit interaction, the majority regime in the phase 
diagram is gapless with pairing symmetry of supercon¬ 
ductivity being p + d-wave. Furthermore, the gap func¬ 
tion d -vector for p-wave superconductivity is not aligned 
with the internal magnetic field of the spin-orbit interac¬ 
tion. As the spin-orbit interaction is turned on, we find 
that the ground state undergoes a phase transition to 
a topological gapless phase with each nodal point origi¬ 
nated from pure d -wave being split into two stable nodal 
points characterized by the symmetry class DHL Due 
to the splitting nodal structure, zero-energy Majorana 


modes always exist for any interfaces that are not ex¬ 
actly in (100) or (010) directions. In addition, due to 
non-aligned d- vector, as the p-wave amplitude further in¬ 
creases beyond the Rashba spin-orbit interaction, a tran¬ 
sition of nodal point configuration occurs. The ground 
state is still gapless. However, the Majorana flat band 
near p y = 0 in (110) edge becomes the dispersive Majo¬ 
rana edge modes. 

In addition to the mean-field phase diagram, the tun¬ 
neling conductances for the NIS and SIS' junctions are 
also computed. Our results indicate that due to the pres¬ 
ence of dispersive Majorana modes, a small plateau with 
shoulders near the zero bias peak would be induced in the 
tunneling spectrum. This small plateau with shoulders 
is a result due to the Dirac cone of dispersive Majorana 
modes at the edge. The appearance of the Dirac cone 
at p y = 0 shown in Fig. 11 is a manifestation of p-wave 
order parameter and results from competition between 
the d-wave and p-wave pairing symmetries. When the d- 
wave dominates for |or| > |A p |, there are only flat bands 
at the edge, exhibiting a single ZBCP in the tunneling 
spectrum. Only when \cxr\ < |A P |, the Dirac cone near 
p y = 0 emerges and results in the additional feature that 
exhibits the feature of a small plateau with shoulders. 

The existence of the Dirac cone at the edge requires 
particular symmetries. General perturbations such as 
disorders do not respect these symmetries and would not 
yield the feature of a small plateau with shoulders in 
the tunneling spectrum 38 . Hence the corresponding fea¬ 
ture can not be generated by disorders. The feature of 
a zero-bias conductance peak sitting on a small plateau 
with shoulders has been frequently ob serve d in tunneling 
spectrum of (110) edge during the pastP^^. Even though 
these evidences may not be the unique signature to con¬ 
clude the existence of dispersive Majorana fermions and 
further evidence from other measurements is in order, 
they show supportive evidence for dispersive Majorana 
modes and indicate the importance of the spin-orbit in¬ 
teraction in high-T c cuprates. 

In addition, we find that the overlap of dispersionless 
Majorana flat-bands in a Josephson junction always re¬ 
sults in 47r-periodicity in Josephson effect. In particular, 
for typical configurations of 7r-ring in tricrystal experi¬ 
ments that can trap a half-flux vortex, we find that as 
long as the difference of the orientations for a axes of two 
crystals involved in a junction is within 21° — 39°, the 
junction always exhibit 47r-periodicity in Josephson ef¬ 
fect. For typical configurations of tricrystal experiments, 
one has 0a = 30°, 6 = 56.5°, and 0b = 56.5°^. The 
edge of crystal B in facing A is along (010) direction. If 
the orientation of interfaces is perfect, since there is no 
edge state in (010) edge, one expects that there is no Ma¬ 
jorana fermion trapped in the 7r-ring. However, in real 
experiments, there may exist deviations in orientations 
of edges. Any deviation of the edge for crystal B from 
the (010) direction would satisfy Eq.(30). This would re¬ 
sult in trapped Majorana fermions between A and B so 
that the trapped flux in the 7r-ring jumps in unit of two 
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flux quanta in externally applied magnetic fields. This 
would be consistent with the experimental observation 
that 3/2-flux quanta has not been observed in tricrystal 
experiments 31 . Therefore, the above analyses indicate 
that Majorana fermions may have been already observed 
in tunneling experiments on high-T c cuprates. 

While so far in this work we only consider results based 
on the mean-field theory, we expect our results are ro¬ 
bust in the presence of correlation effects as long as the 
symmetry of the system is not changed. Even though 
our results show agreement with past experimental ob¬ 
servations, definite confirmation of Majorana fermions 
in high-T c cuprates requires further experimental stud¬ 
ies. Nonetheless, our results offer important signatures 


to look for in future experiments. In particular, the 
setup configuration for tricrystal experiments would of¬ 
fer a unique way to hold Majorana fermions in high-T c 
cuprates. All the signatures we find are crucial for suc¬ 
cessfully searching Majorana fermions and are left for 
future experimental confirmations. 
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